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1 ABSTRACT 

|Ve compute the variation of the beaming fraction with the efficiency of high energy 7-ray 
production in the outer gap pulsar model of Romani and Yadigaroglu. This allows us to 
correct the fluxes observed for pulsars in the EGRET band and to derive a simple estimate 
of the variation of efficiency with age. Integration of this model over the population of young 
neutron stars gives the expected number of 7-ray pulsars along with their distributions in age 
and distance. This model also shows that many of the unidentified EGRET plane sources 
should be pulsars, and predicts the 7-ray fiuxes of known radio pulsars. The contribution of 
unresolved pulsars to the background fiux in the EGRET hand is found to be about 5 %. 
Subject headings: pulsars — gamma rays — statistics 



1. Introduction 

To date, the only identified galactic point sources of 
high energy 7-rays are young, spin-powered pulsars. Five 
pulsars have been now detected at energies E > 100 MeV 
by the iJC-R-E'T experiment on CGRO, several other sources 
have plausible pulsar counterparts, and a number of ad- 
ditional pulsars have been detected at lower energies. Of 
the EGRET sources one, Geminga, is not detected as a 
radio pulsar; conversely several otherwise likely radio ob- 
jects have not been detected at EGRET energies. Further, 
there are ~ 30 — 40 sources that have been detected by 
EGRET at low galactic latitude (Fichtel et al. 1994). It 
is important to estimate the fraction of these sources that 
are rotation-powered pulsars similar to the observed ob- 
jects. 

In the outer magnetosphere model of Romani and Yadi- 
garoglu (1995; hereafter RY), we developed a model for the 
beaming of high energy 7-ray emission that reproduces the 
observed properties of individual 7-ray pulsars, including 
high energy pulse profiles, the relation to the radio pulses, 
and the non-detection of pulsars such as PSR B0656-I-14. 
In RY, it was shown that the statistics of 7-ray pulsars 
and Geminga-like objects were generally in good agree- 
ment with the predictions of our 7-ray pulsar model. Our 
results are also in general agreement with earlier popu- 
lation studies based on extrapolation from the observed 
pulsars (c/. Helfand 1994), but because of our detailed 
model predictions our inferences can be made much more 
specific. In this paper we examine the population predic- 
tions in detail and compare these results with the observed 
galactic plane sources. 
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To this end we derive the evolution of the 7-ray beam- 
ing factor and the pulsar luminosity with age. This evo- 
lution differs substantially from that of other high-energy 
pulsar models and shows that the 7-ray pulsars are a bi- 
ased subset of the young pulsar population. In addition 
these sums allow us to estimate the mean distances and lu- 
minosities of the galactic plane sources. Comparison with 
the observations shows that most of these are likely to 
be young pulsars with lO'* < < 10® y. These computa- 
tions give important guidelines to assist in detecting other 
young pulsars and provide estimates of the contribution 
of these objects to the general fiux in the galactic plane. 
Finally, there are other widely studied models for the pro- 
duction of pulsar 7-rays in the literature. The most im- 
portant of these is the "polar cap" model (c/, Harding and 
Daugherty 1993, Sturner and Dermer 1994). We compute 
the statistics and expected luminosities of 7-ray pulsars 
in this model, finding that in the context of the standard 
picture of radio pulsars these models do not reproduce the 
properties of the observed population, in contrast to our 
outer magnetosphere sums. 

2. Beaming Factor and Efficiency Evolution 

In RY we showed that the 7-ray pulses observed from 
young pulsars can be explained by emission arising from 
the upper surface of an acceleration gap in the outer mag- 
netosphere, when the full effects of aberration, retarded 
potential and time-of-fiight across the magnetosphere are 
included. This model is an outgrowth of the "outer gap" 
picture of Cheng, Ho and Ruderman (1986, CHR) (see 
also Ho, 1989), but entails the full zone extending from 
the null charge (SI • B = 0) surface to the speed of light 
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Fig. 1. — The beaming factor as a function of the to- 
tal gap efficiency rj (including both photons and pairs) 
for various a. The diagonal lines are combinations of 
and rjj allowed by the measured photon flux for individ- 
ual pulsars. The a range indicated by polarization and 
beam shape studies is shown with heavy lines. Geminga's 
distance is in the range 250 to 500 pc; our efficiency law 
gives d ~ 400 pc; the resulting a agrees with the 7-ray 
pulse shape. 

cylinder, bounded below by the surface of last closed field 
lines, and above by field lines a variable distance away. 
The gap width is specified at the null charge surface and 
incorporates some rough idealizations of the physical pro- 
cesses limiting gap growth. 

In CHR it was argued that as a 7-ray pulsar ages it 
will experience greater difficulty in completing the pair 
production that establishes the upper surface of the gap 
and the radiation zone. We model here a similar behav- 
ior and note that as the gap widens and the distance be- 
tween the last closed field lines and the radiating surface 
increases there are two important effects. First, more 
of the open field lines occupy the gap; this implies that 
a larger fraction of the spin-down power is expended in 
accelerating charges at ~ the co-rotation charge density 
Pgj = 0.07 B12 cm~^ in the gap. This power is 

down-converted and re-radiated at the upper surface of 
the gap zone (Chiang and Romani 1994) and thus older 
pulsars with larger gaps will be more efficient at produc- 
ing ~ GeV 7-rays (CHR). Second, as the gap widens, the 
locus of the radiation surface moves closer to the dipole 
axis and the 7-ray beam narrows. This results in a smaller 
beaming fraction with the 7-rays largely concentrated to 
the spin equator. 



The gap thus spans a variable fraction of the open 
field lines; with an assumed constant radiation efficiency 
from the accelerated particles this fraction is in fact rj^ = 
E-y/EsD = Ej/lilil, the 7-ray efficiency. For each gap 
width we compute both this efficiency and the fraction 
of the solid angle swept out on the sky by the result- 
ing 7-ray beam. The true gap efficiency should include 
the energy emitted in the relativistic wind of pairs; for 
the Crab this is E^± = N^±jmeP ^ lO^^-^^s"! • 10^(2 x 
511 keV) ~ 5 X 10^^ erg s"-'^ ^ Lj. So for steep spectrum 
pulsars with copious pair production such as Crab and to 
a lesser extent Vela, the true gap width is sufficient to en- 
sure that the real efficiency rj remains larger than ~ 10 %. 
While this is important to the details of the pulse shape 
for these very young pulsars, it makes little difference to 
/-y . In Figure 1 we plot against rj for a range of mag- 
netic inclinations a. The decrease in beaming factor with 
increased efficiency is clear, although previous works (cf. 
Helfand 1994) have assumed that is uncorrelated with 
rj and in fact constant. 

In RY we showed that the a of the observed 7-ray pul- 
sars can be estimated from both the 7-ray and radio data. 
In addition, the observed EGRET 7-ray fiuxes and dis- 
tances allow us to compute the pulsar luminosity phase- 
averaged wtthm the j-ray beam. Since in this model, 7-ray 
beams have a roughly constant phase-averaged fiux, we see 
that the total 7-ray luminosity is 

= Aird'^E^ ■ = r]^ ■ Esd- 

Thus the observed fiuxes E^, distances and spin-down lu- 
minosities define the diagonal lines shown in Figure 1. 
With the observational estimates of a, we can then de- 
rive the true 7-ray efficiencies. These are given in Table 1. 
While the pulsed luminosity peaks strongly at ~ 1 GeV for 
the older objects. Crab has a much softer spectrum and it 
is important to include the X-ray fiux in its total emitted 
power. We estimate that the fraction of power in the GeV 
band scales as T4 for pulsars with a characteristic age less 
than 10'*r4 y. From the observed values for PSR B1706-44 
and PSR B1055-52 alone we derive the phenomenological 
evolution of rj^ for all high energy pulsed emission, finding 

7]-^ = 3.2 X 10-5 r°^^« (1) 

with characteristic age r in years; the results agree well 
with estimates for the other objects. Note that the beam- 
ing corrections result in a significantly fiatter evolution 
curve than the r~^^ fitted by Thompson ei al. (1995, 
T95). Detailed models of gap closure will be needed to 
replace this empirical law by the theoretical scaling with 
spin parameters. The cut-off to this evolution occurs as 
the beaming factor drops to zero at ~ 10® y. 

We have ranked the known pulsars by estimating their 
E > 100 MeV fiux from (1) and their distances, using an 
average inclination a = 60°. The top 15 objects for which 
/-y is also > 5 % are listed in Table 1; five of these are 
the observed EGRET pulsars. There are EGRET point 



Submitted to the The Astrophysical Journal, 25 October 1994 



3 



Table 1: 7-Ray Pulsars and Candidates: Fluxes, Efficiencies and Beaming Fractions 



PSR 


off Tl 7/ 1 


d(Vuc) 


P a 
rj 


"I ) 


f 

J J 


^ mod 


1 1'y V_y W 111 111 <3 11 Lj 


B0531+21 


3.10 


2.00 


10. 


80 


0.86 


9.8 


0.0072 ~ 110 including X-ray flux 


B0833— 45 


4.06 


0.50 


71. 


65 


0.76 


120. 


0.038 


BIyOd— 44 


4.24 


1.82 


8.3 


45 


0.55 


8.3 


0.053 


J0d30+17 


5.53 


~0.40 


37. 


~40 


0.23 


37. 


0.51 


B1055— 52 


5.73 


1.53 


4.2 


~90 


0.18 


4.2 


0.71 


B1951+32 


5.03 


2.50 


<1.7 




0.53 


20. 


Wide pulse suggests small a 


B0535+28 


5.78 


1.50 






0.15 


10. 


EGRET somce J0542-|-26f -F-y ~ 2.8 


T>ni 1 A 1 

r)Uii4+0o 


K A A 
0.44 


/.i/ 


<i.O 




U.o/ 


0.0 




B0656+14 


5.05 


0.76 


<1.1 


-10 


0.25 


4.8 


l3 = +8° , beamed away 


B0355+54 


5.75 


2.07 


<1.9 


50 


0.16 


4.1 


/3 = +4°, poor a fit, pulsed fiux? (T95) 


B0740-28 


5.20 


1.89 


<1.9 


40 


0.32 


2.9 


l3 = —10°, beamed away 


B1221-63 


5.84 


2.29 


<2.2 


50 


0.10 


2.6 


l3 = +7° , beamed away 


B1634-45 


5.77 


3.83 


<3.9 




0.15 


2.2 




B1823-13 


4.33 


4.12 


<3.4 


-30 


0.42 


2.1 


EGRET souice J1823-12f ~ 9 


B1046-58 


4.31 


2.98 


<2.2 




0.70 


1.6 


EGRET source J1047-58, 7 



^E-j > 100 MeV 7-ray flux in units of 10 ^^( ^'^^ ). For the candidates we convert from EGRET count rate limits on narrow pulses (T95) 

using {E-^) = 750 MeV, as appropriate to pulsars in this age range. Three pulsars have plausible associations with EGRET point sources. 
^Using Eq. 1 and the estimated a. For the candidate pulsars with measured radio polarization sweep we use the a value from Rankin (1993). 

Where none is available we have assumed a = 60°. 
^RY used a ~ 70 of Lyne and Manchester (1988), which happened to give the correct 7-ray pulse width and phase offset for the generic 

T ~ 10^*^ gap used. With the evolved gap, a ~ 90 (Rankin 1993) is convincing as it gives the observed 7-ray pulse and radio interpulse. 
^The EGRET sources consistent with PSRs B0535-t-28 and B1823-13 may be associated with the supernova remnants G 180.0-1.7 (S147) and 

G 18.8-1-0.3 (Kes 67) (Sturner and Dermer, 1995); much of the flux may be due to local cosmic ray enhancement. 



sources possibly associated with three of the candidate 
pulsars. Only pulse upper limits exist for the other ob- 
jects. In a number of cases, we have estimates of a; for 
these pulsars Fmod and have been corrected accord- 
ingly. It is interesting to note that for several of the pul- 
sars with upper limits below our predicted fiux the known 
magnetic geometry indicates that the 7-ray pulse should 
be beamed away from the Earth. From an examination of 
the /-y, this is likely to be true for most the other sources, 
as well. 

3. Population Estimates 

The beaming factor and efficiency evolution derived 
above allow us to integrate over the population of young 
pulsars. We adopt conventional assumptions, namely that 
pulsars are born at a rate R ~ 1/100 y at spin peri- 
ods of ~ 10 ms, that magnetic fields are distributed as 
a Gaussian in log(B), with mean log(B)=12.7 and disper- 
sion (T(log B) = 0.3, and that the magnetic inclination a 
is distributed isotropically at birth. We ignore any expo- 
nential field decay for these young objects. The evolution 
of the radio beam is taken from Biggs (1990): 

VF(P) = 6.2° p-1/2 

while for the pulsar luminosity, in standard units of 
mJy kpc^ at 400 MHz, we follow Lorimer ei al. 1993 



and take 

= 2.8 mJy kpc'^ e^/^ 

where e = min(PP/10~^®, 10^), with a Gaussian spread 
in log(_L4oo) of (Ti = 0.8. These pulsars are distributed in 
the galaxy with a radial law ^!r(R) R dR that peaks at 
the ~ 5 kpc ring, falling as a Gaussian towards the outer 
galaxy 

f exp(- ^l^l~.f/^' ) X 0.743, |i?| < 3.7 kpc 
[exp(-^), |i?|>3.7kpc 

where i?j = 1.5 kpc, i?o=4.8 kpc and the constant serves 
to match the inner and outer galaxy solutions (Johnston 
1994). We also give the pulsars a Gaussian scale height at 
birth of zq ~ 80 pc and a galactic z velocity drawn from 
the 1-D projection of the 2-D Lyne and Lorimer (1994) 
distribution. 

A Monte Carlo integration of the population consists 
of samples at each Galactic longitude /, with draws of 
the pulsar properties a, B, t, zq and V;,. We compute 
the 7-ray and radio luminosities and beaming factors, de- 
termine the fraction of the sky that is swept out by the 
radio and/or 7-ray beams to find the detection probability 
and the maximum detectable distance. Pulsars are then 
populated at various d in direction / according to the nor- 
malized population ^[i?(/, d)] and then placed in the h bin 
corresponding to the height \z\ = \zo + v^, ■ t{T)\ to form 
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Fig. 2. — The expected number of detectable 7-radio pul- 
sars (solid line) and Geminga-like pulsars (dashed line) as 
a function of characteristic age r. EGRET pulsars (open 
circles) and Geminga (triangle) are also indicated. 

the Monte Carlo sample. For the radio pulsars, our goal 
is not to model the observations in detail, so we do not 
model individual surveys; a flux threshold of 5 mjy ap- 
proximates the average survey depth for the bulk of the 
young pulsar population. For the 7-ray detections, we 
adopt a flux threshold of 3 x 10"^*^ erg cm~^ s~^, compa- 
rable to the faintest 5(7 sources in the flrst EGRET galactic 
plane source catalog (Fichtel ei al. 1994). 

Several important features of the populations are ex- 
plained by our sums. In Figure 2 the detectable num- 
bers of objects, per decade in r, expected for EGRET am 
displayed as a function of characteristic age. Noting 
that the total number of radio pulsars expected in the 
EGRET band is 5 in this sum, the agreement with the 
observed sample (open circles) is quite good. The number 
of Geminga-like objects not visible in the radio is 17. Only 
Geminga (triangle) has at present been identifled, but sim- 
ilar objects are believed to be present in the unidentifled 
source list. We assess this possibility by comparing with 
the source catalog in the next section. Very deep radio 
surveys would increase the fraction of radio-detected pul- 
sars to ~ 1/2 in our model. Varying the initial B flelds 
within a reasonable range changes the total number of ob- 
jects by a factor of order two. The beaming evolution 
prevents most pulsars as old as 10® y from being seen as 
7-ray objects. 

The limited 7-ray sensitivity assures that the high en- 
ergy pulsars are a more local sample of the galaxy than 
the radio objects, as is shown in Figure 3. Our sums 



also explain that, like the observed radio pulsars, the de- 
tectable 7-ray pulsars are a biased subset of the parent 
pulsar population (cf. Bailes and Kniffen 1992). In partic- 
ular for a given r the observed sample is heavily weighted 
towards low B (and low P). This is a natural result of 
our beaming and efficiency evolution, and is clear in the 
EGRET pulsars. However, it should be noted that very 
low B flelds (e.g. PSR B1951-F32) may not allow effi- 
cient particle multiplication and GeV 7-ray production, 
while high held objects (e.g. PSR B1509-58) may have a 
strongly enhanced soft synchrotron spectrum; further po- 
larization studies are needed to see if these objects exhibit 
true spectral variations or primarily beaming effects. 

Looking at the flve EGRET pulsars it can be speculated 
that some initial B held decay occurs on short time scales. 
We have run our simulation with a power law B held decay 
that matches the trend observed for the EGRET pulsars. 
The effect is to increase the number of detections for r 
between 10® and 10'^ y. 

Before we leave the population estimates, we note that 
we have computed the population models for the luminos- 
ity and beaming laws of the polar cap model of Dermer 
and Sturner (1994). Using their formulation, we flnd that 
only 0.04 radio pulsars and 0.03 radio quiet objects should 
be detected in the present EGRET source catalog, and less 
than 0.01 objects as young as the Crab are expected. An 
error in Dermer and Sturner's sums was the assumption of 
a source detection threshold of 10""'^"'^ erg cm~^ s""'^; such 
fluxes may be reached at ~ 3(7 at high galactic latitude 
for long integrations, but are much fainter than the source 
detection threshold in the galactic plane. Of course, if we 
assume increased iJC-R-E'T sensitivity, the properties of the 
detected objects match poorly with those of the observed 
pulsars in this model; in particular, comparing with Fig- 
ure 3, the mean source distance of ~ 10 kpc is seen to 
be much larger than that of the observed pulsars. Other 
polar cap models fare similarly poorly in population sums. 

4. Unidentified EGRET Sources and Conclusions 

In the EGRET Galactic plane (|&| < 10°) catalog 
there are 42 sources with detection signiflcance > 5(7, flve 
of which are the young pulsars mentioned above. Our 
goal is to estimate what fraction of this sample is due to 
unidentifled young pulsars. Several caveats are in order. 
Given the extreme sensitivity of the source detection to 
the background model, it is likely that a modest fraction 
of these sources are spurious. Conversely several addi- 
tional steady sources will be found at comparable levels in 
later catalogs. Present in this sample are certainly some 
AGN visible through the plane of the Galaxy; extrapo- 
lation from high latitude counts suggests that this may 
contribute ~ 10±5 sources to the plane sample (Kanbach 
1995). Additionally, concentrations in the diffuse emis- 
sion, either from molecular clouds or from SNR enhance- 
ments in the local cosmic ray density (Sturner and Dermer 
1995) may also be present, if unresolved. Nevertheless, we 
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Fig. 3. — The computed median distance and B field for radio pulsars (dotted line), 7-radio pulsars (solid line) and 
Geminga-like pulsars (dashed line) as a function of characteristic age r. Observed objects are also indicated: young radio 
pulsars (dots), EGRET pulsais (open circles). Table 1 candidate pulsars (small squares), and Geminga (triangle). 



can use the available source positions and fluxes to check 
agreement with our model predictions. 

In Figure 4 we follow Kanbach (1995) in showing the 
distributions in / and h of the EGRET sources (his- 
tograms) along with the model predictions for radio- 
detectable and Geminga-like 7-ray pulsars. We have as- 
sumed a symmetric galaxy and reflected the data about 
/ = and 6 = 0. The agreement confirms the average 
source distance in our sums. We note that there is an 
excess at low h compared to the model; this is likely due 
to contamination by "diffuse emission" sources. We have 
also found a correlation between the mean h angle and age, 
due both to proper motion and to the decreased sample 
depth for these fainter objects. Finally, from our sums we 
also expect one or two unidentified sources to be pulsars 
with l&l > 10°. 

What may be done to further test this model? Fore- 
most, spectral data on the unidentified sources will be 
essential in separating the pulsar candidates (with spec- 
tral indices ~ —2 to — 1.3) from the much softer galac- 
tic emission. To a lesser extent AGN may be identified 
by spectral differences, but these sources should be pri- 
marily separated by 7-ray variability (present in virtu- 
ally all high latitude sources). The remaining hard spec- 
trum, steady sources are good pulsar candidates. A few 
of these objects should yield counterparts to very deep ra- 
dio searches. Also several pulsar IDs should be supported 
by identification of plerionic SNR in the EGRET error 
boxes. To some extent the &(r) correlation expected for 
pulsars may be tested in the unidentified sources by noting 
that the detected pulsars show some spectral index vari- 
ation with r. Duplication of Halpern and Holt's identifi- 



cation of lower-energy (e.g. X-ray) pulsations, will prove 
difficult except for the brightest sources. Scaling from 
the Crab and PSR B1509-58, ASCA observations for the 
youngest pulsars may have some success. The sample of 
both Geminga-like objects and 7-ray/radio pulsars can of 
course be increased by more careful analysis of the obser- 
vations, allowing detections to fainter fiux levels. 

Looking at the candidate pulsars in Table 1, one can see 
that if further analysis of the increasing EGRET photon 
database can decrease the limiting sensitivity by a factor 
of ~ 3 — 5, a number of additional radio pulsars should 
be within reach. Given the dependence of the model fiux 
on a and corrections to the fiux from the details of the 
pulse shape, a number of pulsars ranked below those in 
Table 1 will also be interesting targets; an additional ~ 20 
objects have model fiuxes within a factor of 3 of that for 
PSR B1823-13. However, at typical ages of log(r) > 5.5 
most of these will have small and will not be visible 
from Earth. Since the 7-ray beams are confined to the spin 
equator for the efficient older pulsars; we require ( = a + 
[3 ^ 80° to see these pulsars. Objects with interpulses will 
have large a and are thus excellent targets; PSRs B1719- 
37, B1736-29, B0906-49 and B1822-09 are likely to be seen 
with a factor of 5 increase in EGRET pulse fiux sensitivity. 
We note also that our beaming model provides expected 
7-ray pulse phases and profiles, which can enhance the 
significance in these pulse searches. 

We can at some level test the significance of the beam- 
ing corrections from Table 1. For r ~ lO^y pulsars, /-y 
suggests that ~ 80 % should be detected; 1/1 are seen. 
For T ~ lO'^y pulsars the detection fraction should be 
~ 60 %. At present 2/4 from Table 1 are detected; we 
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Fig. 4. — The distribution of unidentified EGRET galactic plane sources in / and h. Curves show model results for 7-ray 
pulsars with (solid line) and without (dashed line) detectable radio beams. The histograms of identified EGRET sources 
and the observed pulsars are also shown. A total of 10 AGNs distributed uniformly were added to the dashed curves. 
The shaded bin shows the total number of expected pulsars for |&| > 10°. 



might expect one of the other candidates to be confirmed. 
For 10® - lO'^y pulsars, the mean /y is ~ 25 %. Of the 
9 objects in the Table, 2 are presently seen (not counting 
Geminga); we might therefore expect one of the candi- 
dates to be confirmed. We have not included an a bias 
implicit for radio detections in these estimates. 

Finally, it is traditional to estimate the contribution 
of the pulsars not visible as discrete sources to the dif- 
fuse galactic background (e.g. Bailes and Kniffen 1992). 
In our model, the integrated fiux from sources below our 
identification threshold is 1.9 x 10~® erg cm~^ s~^. This 
is only about 5 % of the measured plane fiux. To some 
extent this may be testable since this component should 
have the E~^^^ spectrum and / and h distribution of the 
T ~ lO^y pulsars. 

In conclusion, it is clear that the outer magnetosphere 
model of RY can provide a good explanation for much 
of the galactic 7-ray source population. While the lumi- 
nosity law is at present phenomenological, the rj^ and 
variations in the model are important in interpreting the 
observations. Further observations and analyses can test 
this model and more careful studies of the brightest pul- 
sars, coupled with detailed spectral modeling (Chiang and 
Romani 1994) can illuminate the mechanics of the pulsar 
magnetosphere itself. 
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